1/ f a spectra in elementary cellular automata and fractal signals 
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We systematically compute the power spectra of the one-dimensional elementary cellular automata 
introduced by Wolfram. On the one hand our analysis reveals that one automaton displays 1/f 
spectra though considered as trivial, and on the other hand that various automata classified as 
chaotic/complex display no 1/f spectra. We model the results generalizing the recently investigated 
Sierpinski signal to a class of fractal signals that are tailored to produce l/f a spectra. From 
the widespread occurrence of (elementary) cellular automata patterns in chemistry, physics and 
computer sciences, there are various candidates to show spectra similar to our results. 



PACS numbers: 05.45.Df, 89.75.Da, 82.40.Np, 45.70.Qj 

In 1984 Wolfram introduced the so-called elementary 
cellular automata (EC A), opening a field still being 
vividly active 20 years thereafter jjj. Wolfram's more 
recent popular book Q has attracted great attention, 
although the opinion of the work's merits are divided 
among the scientific community ,3j. EC A are discussed 
extensively in the context of computationally irreducibil- 
ity of physical systems |4|, e.g. it is proven that in the 
Turing sense [5| rule 110 (being one of the possible 256 
EC A) is an universal computer Moreover, possible 
transformations between difference equations and (E)CA 
have been investigated . Among the numerous physical 
applications we mention here only (kinetic phase transi- 
tions in) catalytic reaction-diffusion systems 0, 0, El i 
deterministic surface growth 11], branching and annihi- 
lating random walks 12] and random boolean networks 

It is important to note that Wolfram's ECA are often 
studied for different boundary conditions on a finite ar- 
ray. A particular boundary condition (e.g. a periodic 
or an absorbing one) disturbs the pure evolution of an 
ECA. As a result, some automata display complex be- 
havior, while other are simply periodic. Though there is 
no algorithm for classifying a given elementary automa- 
ton, Wolfram conjectured that ECA can be grouped into 
four classes of complexity: 

Class 1: Steady state, class 2: Periodic or nested struc- 
tures, class 3: Random (chaotic) behavior, class 4: Mix- 
ture of random and periodic behavior. 

The first class represents automata that are (for al- 
most all initial conditions) trivial in the sense being static 
or finally evolve to the some steady state. Those rules 
that belong to the second class produce simple periodic 
or self-similar, i.e. fractal, structures. The third class in- 
cludes rules exhibiting random patterns, e.g. a particular 
rule (number 30) is used to generate random numbers in 
Mathematica. The fourth class is somehow a mixture of 
classes 2 and 3 generating the most complex structures. 
For more rigorous classifications we refer the reader to 



the literature [lll^]. 

Since the coining paper of Bak, Tang, and Wiesenfeld 
|l5| , there has been considerable interest in the long-time 
behavior of cellular automata, especially for occurrence 
of long range correlations, and correspondingly for power 
spectra exhibiting a power law decay S(f) ~ /~ Q with 
a w 1.0. Despite the abundance in nature, systems ex- 
hibiting spectra with exponents near to 1 are poorly un- 
derstood. While the mechanisms generating 1/ f a spec- 
tra may be substantially different from each other, some 
models and the observed l// a power laws have become a 
paradigm for complex dynamical systems in general [l^ ] 
(see also references in |lCl|). 

Definition of ECA. An elementary cellular automa- 
ton consists of an infinite one-dimensional lattice of cells 
being either black (1) or white (0), and a deterministic 
update rule. At each discrete time step, a cell is updated 
x*j — > a4 +1 according to the state of the next-neighbor 
sites and its own state one time step before: 



x n f i x n+l> X n> X n-l) 



(1) 



where / (the rule) is determined by 8 bits being the out- 
put of the possible input bits 000, 001, 111. As a 
consequence, there are 256 (ECA-)rules that are named 
rule - 255. In this article we focus on rules 90 and 150 
defined by 



b n+l 



mod 2 



(2) 



where r = defines rule 90, and r = 1 rule 150, re- 
spectively. As demonstrated earlier, rule 90 can be in- 
terpreted in the context of catalytic processes. A pro- 
cess (catalysis) is initiated (or continued) when exactly 
1 neighbor site is active whereas the process (catalysis) 
is stopped when too mam/, i.e. 2, or to less, i.e. no, 
neighbor sites are active [lOj. 

A similar interpretation may be given for rule 150. 
Catalysis at x l n is stopped when no or two neighbor sites 
(now x l n included) are active and it is initiated (or contin- 
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ued) when one or three sites are active. Note that both 
rules mimic local self- limiting reaction processes 0, 0] . 

Spectra of sum signals. It is known that ECA on fi- 
nite lattices for various boundary conditions display no 
1// Q spectra Q. Rather than evaluating the rules on 
finite lattices we calculate the evolution on an infinite 
lattice. More precisely we focus on a sum signal defined 
as the total (in)activity, magnetization, etc. of the whole 
system: 



(3) 



We have systematically investigated all 256 rules, for lo- 
calized initial conditions (i.e. single 1, 11, 101, 111, ...), 
as follows. The sum signal for non-trivial rules exhibits 
increasing mean (X) t well fitted by a power law in time 
[2fij | . Consequently, we focus on the detrended sum signal 
defined by 



Y(t) = X{t) - f{t) 



(4) 



where the coefficients of f(t) = at b are fitted. However 
for some ECA Y(t) possesses an increasing mean vari- 
ance. Thus we investigate for each automaton another 
signal (and its spectrum) 



Z(t) = Y(t)/{Y)f_ 



i+i,t+i} 



(5) 



where 21 is the width of a sliding window that normalizes 
the fluctuations of the detrended signal Y{t) according to 
the method of detrended fluctuation analysis (DFA) ap- 
plied for non-equilibrium processes 0]. We have calcu- 
lated the corresponding power spectra |A(w)| 2 , |i^(w)| 2 
and \Z(uj)\ 2 for all 256 ECA. It turns out that that 25 
of the 256 rules exhibit 1//" spectra whereas 231 do not 
(see table HJ. 23 of those automata that exhibit l/f a 
spectra display Sierpinski patterns, i.e. well studied self- 
similar structures |lfj. Their spectra are extensively in- 
vestigated in 0] exhibiting 1// Q spectra with exponents 
1.15 ±0.05. 

The two other rules, i.e. 105 and 150, show a different 
behavior. Here we focus on rule 150 21 1. The first 128 
time steps of the evolution for a single 1 as initial condi- 
tion is depicted in Fig. 2] (upper inset). It is a Sierpinski- 
like self-similar structure. However the fractal dimension 
differs from the Sierpinski gasket (d = 1.58) being the 
golden mean d = (1 + v / 5)/2 w 1.69. Fig. □ shows also 
the corresponding signals X(t) and Z(t). The spectrum 
Y(lu) is displayed in Fig. [21 For u) not too small, the 
averaged spectrum exhibits a straight line in the log-log- 
plot verifying a power law behavior. Depending on the 
average process and fit range we obtain a fit exponent of 
about a — 1.27 ± 0.05. Due to dominating randomness, 
members of classes 3 and 4 typically produce thermal 
l// 2 spectra (see Fig. 



Class 


ECA rule number 


1 


218 


2 


(26, 82, 167, 181), (154, 210) 
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(18, 183), (22, 151), (60, 102, 153, 195), (90, 165), 
(122, 161), (126, 129), (146, 182), 105, 150 
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TABLE I: Rules that produce l/f a spectra. Rules in brackets 
belong to one equivalence class. Rules 105 and 150 (bold) 
produce spectra with power law exponents about a — 1.3. 
All other listed rules exhibit spectra with exponents about 
a — 1.2. The 231 rules not listed are not capable to produce 
l/f a spectra, e.g. most of the spectra display no power law 
decay, or exhibit thermal l// 2 spectra (see Fig. |5J . 
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FIG. 1: The first 128 time steps of the time signal X(t) 
generated by rule 150. Upper Inset: Self-similar structure 
generated by rule 150 for the first 64 time steps. Lower Inset: 
Normalized signal Z(t); the straight line indicates Z — 0. 
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FIG. 2: Rule 150 and rule 110: Averaged power spectrum 
of Y(t) up to T/8 for T = 2 18 using (incommensurable) l.l fc - 
bins, i.e. the fc-th interval is defined by [[1.1*1, [Ll fe+1 l] 
where the brackets [ ] denote upwards rounded integer values 
(ceiling function). The inset shows the rule 150 spectrum, 
averaged using 2 fe -bins, i.e. the fc-th interval is defined by 



[2 fc ,2 



fc+i 



1]. Both averages correspond to a constant 5ui/ui 



ratio. The graphs are well fitted by a power law with exponent 
a — 1.27. The thermal l// 2 decay of rule 110 (grey) as a 
typical member of Class 4 is shown for comparison. 
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Fractal signals produce 1/ f a spectra. All ECA that 
are capable to produce a self-similar structure exhibit 
l/f a spectra. Hence one may naively expect that ev- 
ery (self-similar) fractal structure produces 1/ f a spectra. 
However it is important to note that this is not the case. 
There are many fractals like the Koch snow flake, Can- 
tor dust etc. exhibiting no 1// Q spectra because their 
respective sum signals simply grow exponentially [lflj . 

Rather than a geometric approach we focus on fractal 
signals itself. Thus we no w g eneralize the recently inves- 
tigated Sierpinski signal lOj. As we will see, the gener- 
alized signal is capable to model 1// Q spectra producing 
spectra with continuously tunable power law exponents. 
More precisely, we consider the signal 

X 5 (t) = 2 5 ^ aAt} (6) 

where aj is the jth bit of the binary decomposition of the 
discrete time t=0,l,2,.... For 8 = 1 we have shown re- 
cently both numerically and analytically that the signal 
exhibits 1//" spectra with a close to unity. The special 
ansatz, eq. (jBJ), represents a straightforward generaliza- 
tion of the closed form for the sum signal of the Sierpinski 
pattern produced by rule 90 0] . In the next paragraph 
we show that for deviations from S — 1 the signal can 
produce 1/ f a spectra within a wide range of exponents 
a. 

In analogy to the calculation in Ref . j we calculate 
the periodogram X(oj) of the time signal © analytically: 



E e " 

t=0 



X s (t) 



E 



exp(iw V o j 2 j ) X s ( V (Tj V ) 



{rjo,...,aN — x} 



N-1 



e n exp ( iw2j + 6 in 2 )) 

{tr ,...,trs-i} 3=0 
N-1 

= n e exp ( iw2j + 5 in 2 )) 

N-1 

= Yl C 1 + cxp(iw2 J ' + 8 In 2)) . (7) 

3=0 

The absolute value of X(u>) simplifies to a trigonometric 
product which the logarithm converts into a sum: 

iV-l 

In \X{uj)\ 2 = ln[l + 2 2<5 + 2 1+s cos(a;2»)]. (8) 

3=0 

We roughly estimate the sum in eq. JSJ) replacing the sum 
by an integral, and substituting y — w2- ? , 

ln\X(uj)\ 2 » / ln[l + 2 25 + 2 1+,5 cos(w2 i )]dj(9) 



ln[l + 2 2a + 2 1+d cosy] 
^hi2 



dy.(10) 



As ln(a + 6a;) w ln(a) + -x for \x\ <C 1, we obtain 
ln(l + 2 25 ) r 2N ^dy 

^2 L 

cos(y) 



ln|X(w) 



In 2 J u y 

2l+S rU**' 



(l + 2 25 )ln2 



y 



-Ay. 



(11) 



The integral over the integral cosine is nearly indepen- 
dent of the upper boundary for high values of the bound- 
ary. Thus, we can substitute the upper boundary ui2 N ~ 1 
by some A-dependent constant, say cn 3> 1. Finally, 
replacing the cosine by one yields immediately a rough 
approximation of the power spectrum: 



\X(U3)\'< 



C N 0J (1+2") In 2. 



(12) 



For a given power law exponent 0<a<l/ln2« 1.44, 
we obtain 8 from eq. (|12J) as 



In 



8 = 



l + ^/l-Q 2 (ln2) i 
ex In 2 



In 2 



(13) 



To generate signals with goal exponents, e.g., ot\ — 
0.8, a 2 = 1.0, a 3 = 1.2, one can use the correspond- 
ing value of 8 according to Eq. (|13fl . Fig. [3] shows the 
spectrum of the detrended signal © for 62 = 1.31184 
(corresponding to oti = 1.0). For 8\ = 1.72425 and 
8s = 0.902749 the individual spectra exhibit similar 
graphs (not shown). Depending on the averaging the 
power law fits giving ol\ — 0.8 ± 0.1, a\ = 0.95 ± 0.05 
and a3 = 1.2 ± 0.05, are in good agreement with the 
theoretical results. 
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FIG. 3: Averaged power spectrum of the detrended signal © 
for 8 = 1.31184 up to T/8, T = 2 20 using (incommensurable) 
l.l fc -bins, i.e. the fc-th interval is defined by [[1.1*1, [l-l fe+1 l] 
where the brackets [ ] denote upwards rounded integer values 
(ceiling function). The inset shows the spectrum, averaged 
using 2 fe -bins, i.e. the fc-th interval is defined by [2 fc ,2 fe+1 - 
1]. Both correspond to a constant 8ui/io ratio possessing fit 
exponents of about a — 0.93. The corresponding spectrum for 
the variance detrended signal exhibits power spectra around 
a — 1.0. The theoretical exponent is a = 1.0. 
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Two-dimensional automaton. While one-dimensional 
experimental setups as in 01 seem to be quite artificial 
for (self-limiting) catalytic processes, two-dimensional 
dynamics is more generic. Consider the Sierpinski dy- 
namics on a (i, j)-plane: 



;i _ x ] mod 2. (14) 



For a single 1 as initial condition on a plane the sum 
signal X 2 £>(t) — J2i x ij generates the sequence 



X 2D (t) = 1,4,4,16,4,16,16,64,. 



(15) 



More precisely, the recurrence relation generating eq. 
(0 is given by u n — > u n+1 = (u n ,4u n ) for u = (1). 

First, if the factor 4 is replaced by 2, the relation be- 
comes equivalent to the ld-Sierpinski signal A"n>(f) in 
Ref. 0. Second, we obtain X 2 D(t) — Xirj(t) 2 and 
therefore X 2 £,(t) — Xs =2 {t). Thus the generalized Sier- 
pinski pattern in two dimensions exhibits 1/ f a spectra 
with exponents around the value according to eq. i|12|) 
for 5 = 2, that is a = 0.679. We numerically verified the 
value obtaining exponents around a — 0.7 as expected. 

Conclusions. Elementary Cellular automata are a 
paradigm for emergence of complex spatiotemporal be- 
havior from extremely simple dynamics. We systemati- 
cally investigated all 256 elementary cellular automata. 
As expected, among those as (nested) periodic/chaotic 
classified rules (classes 2 and 3) there are various rules 
that display l/f a spectra (see table HJ. Unexpectedly, 
on the one hand all rules classified as complex display no 
l/f a spectra, while on the other hand, the trivial rule 
218 does (being a member of class 1). It is important 
to note that the numerically calculated spectra are ro- 
bust against noise, that is, the fit exponents change only 
slightly for other initial conditions than a single seed. 

Moreover we generalized the approach of a sum signal 
introduced in [lOj to derive analytically the spectra of 
the 2D Sierpinski automaton. The investigated fractal 
signals © serve also as a fit model for signals produced 
by elementary cellular automata rules. We have obtained 
a time series generator with continuously tunable power 
law decay exponent. The tailored signals represent ana- 
lytically tractable (nontrivial) 1//" generators that shed 
light on the arcane mechanisms of 1//" spectra. 

From our results, we expect that in experimental sys- 
tems showing spatiotemporal pattern formation similar 



to the ECA patterns, the power spectra of the total 
(in) activity will exhibit power law behavior within a cer- 
tain range. 
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In [l(j we have shown this both numerically and analyt- 
ically for rule 90. For other rules it is also easy to derive 
analytically. 

Rule 105 is simply the inverse of rule 150, i.e. 
fi 05 (a,b,c) = 1 - f 150 (a,b, c). 



